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Abstract 

With the aim of setting constraints for the modeling of the QCD phase diagram, the phase structure 
of the two-flavor Polyakov-loop extended Nambu and Jona-Lasinio (PNJL) model is investigated in the 
range of imaginary chemical potentials (m) and compared with available Nf = 2 lattice QCD results. The 
calculations are performed using the advanced nonlocal version of the PNJL model with the inclusion of 
vector-type quasiparticle interactions between quarks, and with wave-function-renormalization corrections. 
It is demonstrated that the nonlocal PNJL model reproduces important features of QCD at finite fj,i, such as 
the Roberge-Weiss (RW) periodicity and the RW transition. Chiral and deconfinement transition tempera- 
tures for Nf = 2 turn out to coincide both at zero chemical potential and at finite m. Detailed studies are 
performed concerning the RW endpoint and its neighborhood where a first-order transition occurs. 
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I. INTRODUCTION 



One of the fundamental topics in the physics of the strong interaction is the investigation of the 
QCD phase diagram. Lattice gauge theory provides a powerful tool for dealing with the thermo- 
dynamics of QCD. However, lattice QCD (LQCD) simulations are not feasible at real values (/xr) 
of the quark chemical potential(s), at least not for /xr/T > 1 at given temperature T, because of 
the well-known sign problem. Explorations into the domain of large /x R , or high baryon densities, 
can so far only be pursued in terms of models, a frequently used one being the Polyakov-loop- 
extended Nambu and Jona-Lasinio (PNJL) model yJ-lSD- Such approaches have obvious intrinsic 
limitations. (For example, important constraints well known from nuclear physics around and 
above normal nuclear matter densities are not (yet) realistically incorporated in such models). 
Their predictive power concerning possible phase transitions at high baryon densities is quite lim- 
ited and involves substantial ambiguities at finite /i R [5, (J. Setting additional constraints for the 
modeling of QCD phases at nonzero chemical potentials is therefore an important issue. 

Under such circumstances, the continuation from real to imaginary chemical potential turns out 
to be a useful conceptual strategy. In this case LQCD simulations do not face the sign problem; 
there is no principal restriction in performing LQCD computations in the region of imaginary 
chemical potentials. Models formally designed for general complex chemical potentials, 

can then interpolate between positive and negative /x 2 and incorporate LQCD constraints set at 
imaginary /x into extrapolations to the physical region of real /x. 

The relationship between /xr and /xi is apparent in the connection between the grand canonical 
partition function, Z GC (T, 6), written for imaginary chemical potential in terms of the dimension- 
less ratio 6 = \i\jT , and the canonical partition function Zq(T, N q ) expressed as a function of 
quark number N q : 

Z C (T, N q ) = r ^- e - iN « d Z GC (T, 9) . (LI) 

J-7T 271_ 

Inverting this Fourier transform gives the grand partition function at finite /x R through the fugacity 
expansion as 

oo 

Z gc (T,/xr)= e Nm/T Z c (T,N), (1.2) 

N=-oo 
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where the thermodynamical limit is considered, with the three-dimensional volume set to V —> oo. 
Several important features are noteworthy in the region (see Refs. [Q, kJ- J 



14fl for details). 



First, the QCD partition function has a 2n/3 periodicity along the 9 axis, the Roberge-Weiss (RW) 
periodicity II 1411 . Another characteristic property is the RW transition: thermodynamical quantities 
and order parameters have nonanalytic behavior at 9 = nk/3 (with k an odd integer) once the 
temperature exceeds a specific value, T > Te- This point, (T, 9) = (Te, irk/ 3), is the so-called 
RW endpoint. These nonanalyticities are induced by the transition between different minima of the 
thermodynamical potential (see Fig. 2 in Ref. 111411 ). Moreover, the imaginary chemical potential 
is gauged into the temporal boundary condition of quarks; the dual quark condensate Ill5l - l20ll 
is then proposed as a measure for exploring the possible correlation between the chiral and the 
deconfinement transitions. This correlation was also studied previously in terms of a Ginzburg- 
Landau analysis [21]. 

In recent years several LQCD investigations at imaginary chemical potentials have been per- 
formed in order to draw conclusions about the ur region, either through analytical continua- 



28[|. These studies were so far restricted to a 



tion [22-26] or by the canonical approach [27, 
small range of /xr. To overcome such limitations, the imaginary chemical potential matching ap- 
proach has been proposed in Ref. II 1311 . The basic aim of this approach is the modeling of the 
QCD phase diagram starting from a suitable effective Lagrangian constrained by direct compari- 
son with LQCD results in the region of imaginary chemical potentials /xi. The implementation of 
the correct QCD behavior in the /xi region is then a necessary (though not sufficient) requirement 
for proceeding to physical (real) chemical potentials. 

The analytic continuation of LQCD data from \i\ to /ir commonly assumes a certain form (e. 
g. polynomial) for the interpolating function. This implies a limiting convergence radius. (See, 
for example, Ref. II29L 13011 for a detailed assessment). The strategy pursued in the present work is 
different. Given an effective Lagrangian starting from QCD and its symmetries, a set of possible 
types of effective interactions between quarks is introduced. To the extent that these interactions 
are not entirely determined just by the physical low-mass meson spectrum, lattice QCD results at 
finite \i\ are useful in providing additional constraints. Once the relevant effective interactions are 
fixed, the continuation to finite /x R is not subject to convergence radius limitations any more. 

The PNJL model is of considerable interest in this context. In previous work QVU it was shown 
that the PNJL model can reproduce the RW periodicity together with other important QCD prop- 
erties such as spontaneous chiral symmetry breaking. The RW periodicity is described by the 



extended Z 3 symmetry defined as 



(1.3) 



where n is any integer; $ and are the Polyakov loop and its conjugate, respectively. In the PNJL 
model this symmetry is realized by construction. It has in fact been shown that the local version 
of the PNJL model is capable of reproducing the LQCD data at finite fj,i, but at the expense of 
introducing scalar-type eight-quark and vector-type four-quark interactions [11], or a Polyakov- 
loop dependent NJL coupling strength 113 ill , in addition to the standard (scalar plus pseudoscalar) 
chiral contact interactions between quarks. 

The local PNJL model still works with a schematic ad-hoc momentum-space cutoff that has no 
foundation in QCD. More direct contacts with QCD have recently been established by introducing 



the nonlocal PNJL model 



32, 



330 and by pointing out its formal derivation from QCD [34 



Nonlocal interactions were also investigated, for example, within the framework of NJL 05I1 
and instanton models Ir36ll . It is instructive to extend the nonlocal PNJL approach to imaginary 
chemical potentials in order to examine its properties in direct comparison with LQCD results. 
Very recently, independent related calculations have been performed [37] in parallel to the present 
investigations, leading to qualitatively similar conclusions. The additional new element in the 
present work is the detailed study of the role of (nonderivative) vector interactions between quarks 
as they emerge from their basic color-current couplings. 

In the present work we restrict ourselves to the two-flavor case. Section II briefly summarizes 
the nonlocal Nf = 2 PNJL model in its latest version 13311 with the inclusion of quark wave- 
function-renormalization effects. Numerical results are presented and discussed in Sec. III. The 
paper closes with a summary and conclusions in Sec. IV. 



H. FRAMEWORK AND FORMALISM 

This section proceeds in several steps, starting with a brief exhibition of the nonlocal PNJL 
model and its thermodynamics, followed by some symmetry considerations. Finally, a possible 
additional vector-current interaction between quarks will be incorporated. Such a Lorentz-vector 
interaction is usually not part of the standard PNJL model. It is expected to play a pronounced 
role, however, in studies extended to imaginary chemical potentials. 
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A. Nonlocal PNJL model 



The generic Euclidean action of the two-flavor PNJL model is 



S 



dr / d 3 x[q(x)(i0 - mo)q(x) + C int ) - 0VU($[A),$[A);T) , 



(11.1) 



where q(x) = (u(x), d(x)) T is the two-flavor quark field, m denotes the current quark mass taken 
in the isospin limit (m = m u = ma), and D u = d u + iA u = d u + i5q A°' a \ a /2 is the color gauge- 
co variant derivative, with SU(3) C Gell-Mann matrices X a . The gauge coupling g is understood to 
be absorbed in the definition of A°' a . The last term in Eq. (|II.1I) is the Polyakov-loop-effective 
potential U, multiplied by volume V and inverse temperature (3 = T~ l , and to be specified later. 
As in previous studies we treat the temporal gauge field A as a constant Euclidean background 
field in the form A A = iA° = A\ A 3 /2 + A\ A 8 /2. Further details are given in Sec. HI CI For 



extensions beyond the mean-field treatment of At, see Ref. 



3911 . The nonlocal generalization 



of the PNJL model is characterized by an interaction Lagrangian featuring nonlocal quark currents 



and densities, as follows [SilS, 



32, 



33 



37D: 



C mt (x) = G 



j*(x)j*(x) + J{x)J(x) 



d 4 z C(z) q(x + z/2) T a q(x - z/2) 



J(x) = / d 4 z F(z) q(x + z/2) -!- q(x - z/2) . 

2k 



(H.2) 
(H.3) 

(II.4) 



The chiral (scalar and pseudoscalar) densities j^x) with a = 0, 1, 2, 3 involve the operators -T a = 
(l,i7 5 r). The overall coupling strength G of dimension [length] 2 is chosen sufficiently large so 
that spontaneous chiral symmetry breaking and pions as Goldstone bosons emerge properly. The 
J(x) introduces additional vector- type derivative couplings with q(x') d /1 q(x) := q(x')(d ll q)(x) — 
(d fl q)(x l ) q(x) together with a scale k so that the effective strength of this term in £ int is G/k 2 . In 
the following we refer to the interaction induced by J(x) simply as a derivative coupling in order 
to avoid confusion with the nonderivative vector interaction that will be introduced in Sec. Ill Fl 

Associated with the currents or densities (|II.3I) and (|II.4I) are nonlocality distributions C(z) and 
F(z). As described in detail in Refs. [a la, 



32, 



33f| . these distributions govern the momentum 



dependences of the mass function and of the renormalization factor that appear in the quark quasi- 
particle propagator, Z(p 2 )(^y ■ p — M(p 2 ))~ l . The Fourier transform C(p 2 ) of C(z) is related to the 



quasiparticle mass function M(p 2 ) determined by the self-consistent gap equation, 



M(p 2 ) = Z(p 2 ) m + aC(p 



(11.5) 



where a is the scalar mean field basically representing the chiral (quark) condensate (qq). The 
Fourier transform J-"(p 2 ) of J-(z) is in turn related to the Z factor representing quark wave-function 
renormalization, 



Z{P 2 



i- -Hp 2 

K 



(H.6) 



where v is the mean field induced by J(x) Il40h . 



B. Nonlocality distributions, quark mass function, and quasiparticle renormalization factor 

The following momentum- space forms of the distribution functions appearing in Eqs. (|II.3I) and 
(111.41) are used in the present work: 

C(p 2 )= / d 4 z exp(-ip -z)C{z)= I (II.7) 
J (p 2 >\ 2 ), 

F(p 2 ) = J d A z exp(-ip • z) 7{z) = exp(^-p 2 4/2) . (Q.8) 

The running QCD coupling a s (p 2 ) determines the asymptotic form of C(p 2 ) while its infrared 
behavior is given a Gaussian parametrization with a characteristic length scale dc- The matching of 
these high- and low-momentum representations at an intermediate scale A determines the constant 
J\f. The distribution ^(p 2 ) has its own length scale (If over which the nonlocality unfolds. As in 
Ref. 113 311 we use dc ~ 0.4 fm and dp ~ 0.3 fm. 

With this input, the resulting mass function M(p 2 ) and the renormalization factor Z(p 2 ) closely 
resemble LQCD results, as shown in Fig.[TJ For the mass function we are guided by the lattice data 



of Ref. 114 lh extrapolated to the chiral limit (m — > 0). Both the mass function and the Z factor 



are gauge dependent. Using Landau and Laplacian gauges in comparison, it is demonstrated in 



|4JJ] that the mass function shows very little variation between these two gauge fixings, whereas 
the gauge dependence of the Z factor is about 20 % in the infrared region. This gauge dependence 
of Z is not crucial here since the wave-function-renormalization effects will turn out to be quite 
small in the present context. In practice we use the Landau-gauge results for orientation. 
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Fig. 1. The p dependence of the quark mass function M(p 2 ) and of the quasiparticle renormalization factor 
Z(p 2 ) resulting from the distributions (III.7I ) and (III- 8b (solid lines). Solid triangles and open circles are 
LQCD data generated with a large quark mass, while solid squares show the extrapolation of the mass- 



function data to the chiral limit (from Ref. 114 ill ). Open triangles are LQCD results for Z(p 2 ) taken from 



Ref. [42]. 



At finite temperature T the squared (Euclidean) four-momentum in C(p 2 ) and J-'ip 2 ) becomes 
p 2 = u 2 +p 2 with fermionic Matsubara frequencies co n = (2n+l)7rT. When Polyakov-loop fields 
and nonzero chemical potentials are included, co n is shifted to cu n — ifi plus linear combinations of 



A 



In the PNJL model the pion mass and its decay constant are used to fit parameters in the NJL 



sector of the Lagrangian. These parameters are taken from Refs. [8, 
and with the inclusion of Z-factor effects, respectively. 



3311 for the case studies without 



C. Polyakov-loop-effective potential 



The Polyakov-loop-effective potential U is used in the form given in Ref. [|8|]: 



rp4 



\^b 2 (T)<P<P + b A (T) ln[l -600 + 4(0 3 + 3 ) -3(<P<Pf 



(11.9) 



where <P and <P are represented as 



* = i 

3 



Al + Aj \ ( Al-Al 

exp l i — — + exp — % 



2T 



2T 



exp % 



A* 



$ = $*. (11.10) 
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The coefficient functions b 2 (T) and 64 (T) are parametrized to reproduce pure-gauge LQCD results 



as described in Refs. |8|, 



3311 . The temperature scale T appearing in 6 2 (T) and 64 (T) is set to 



270 MeV, the critical temperature for the first-order confinement-deconfinement transition from 
LQCD in the pure-gauge limit. Variations of this scale in the presence of dynamical quark flavors 
are of potential importance but will not be considered here. 



D. Thermodynamics 



The mean-field thermodynamical potential Q of the nonlocal PNJL model, including quark 
wave-function-renormalization corrections, is constructed using the Nambu-Gor'kov formalism 
(see Ref. [44] and references therein). The final form is 



(i1.1i) 



where 



d 3 p 
(27r)3 



In 



+ 



a 2 + v 2 



AG 



(11.12) 



j = ±,0 71= — OO ' 

Here a and v are the mean fields associated with the scalar density j of Eq. (111.31) and the 
derivative-vector current J of Eq. (111.41) . respectively, upon bosonization. The first term on the 
right-hand- side of Eq. (|II.12I) involves the quark quasiparticle energies 



(11.13) 



with dynamically generated masses, Mj(p 2 ) = M(p 2 = to 2 i + p 2 ), determined self-consistently 
at each shifted Matsubara frequency with i E {±, 0}: 

A 3 A 8 

w„,± = w„ - %\l ± 



2 2^ 

73' 



(11.14) 



where 8 are the gauge fields forming the Polyakov loop already given in Eq. (III. 101) . Likewise, 
the Z factors are understood as Zi(p 2 ) = Z(p 2 = co 2 i +p 2 ). More explicitly: 



Mi(p 2 ) = Zi(p 2 ) m + aC{p 2 = u^+p 2 ) 



Z,{p 2 



l--F{p 2 = u 2 n/l + p- 



-1 



(11.15) 
(11.16) 
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At finite temperature T the Lorentz invariance is broken by the thermal medium and the inverse 
quark quasiparticle propagator becomes S" 1 (p) = Ao(p) 7oP° + A(p)%p l + B(p) with Aq ^ 
A. Here we assume for simplicity that the difference between Aq and A is sufficiently small 
so that it can be neglected, given that the overall influence of wave-function renormalization on 
thermodynamical quantities is not very significant. 

The subtracted ln(u;^ { + p 2 + ml) piece in Eq. (|II.12I) improves convergence in the summation 



over the Matsubara frequencies. The term i?f ree in Eq. (|II.17I) 

d 3 p 



free 



-AT j -^[ln(l + 3(0 + ^ e -^V 



+ ln(l + 3(«P + <p e -^+(p)) e -^+(p) + e -We+{p) j 



(11.17) 



(with (3 = 1/T) is then introduced for consistency. In these "free" parts the quark energies, shifted 
by the chemical potential /i, are taken with the current quark mass m : 



e±(p) 



P' 



The mean fields o and v are determined by the conditions 

dQ _ dQ _ 
do dv 

Similarly, the Polyakov-loop-background gauge fields are determined through 

dQ dQ 



dAl dA s A 



. 



(11.18) 



(11.19) 



(11.20) 



Fluctuations beyond the mean field - not considered in the present work - are of course im- 
portant but do not change the basic phase transition pattern in a qualitative way. Studies of the 
role of selected fluctuations in the PNJL model have been performed previously |38, 39^ . See also 
Refs. flfl] for the treatment of fluctuations in the Polyakov-loop extended quark-meson model. 



E. Symmetry considerations 



48D)as 



The QCD partition function is known to be a function of fi 2 (see for examples Refs. H22I . 
a consequence of time reversal or CP symmetry which implies invariance under the transformation 
fi — > —fi. Our nonlocal PNJL model maintains this property. For later convenience, we introduce 
a modified Polyakov loop and its conjugate, & and as 

& = e u$ t q, = e -&$ . (H.21) 
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These quantities are invariant under extended Z 3 transformations and hence they are RW periodic. 
Properties with respect to charge-conjugation (C) symmetry are of basic interest in this context. 



As Roberge and Weiss 111411 have shown, the RW periodicity of the QCD partition function implies 
that the C symmetry is also intact at 9 = n/3 modulo 2%/3. The #-even thermodynamical potential 
transforms under C as 

Q{9) -> Q{-6) . (11.22) 

Therefore the C symmetry holds at 9 — (2n — l)7r/3 with integer n. On the other hand, quantities 
that are 8-odd transform under C as 

0{9) -> -0{-9) . (11.23) 

Spontaneous breaking of C symmetry is indicated when O has a nonzero expectation value at 
9 = (2n — l)n/3, associated with Z 2 -symmetry breaking under yU-reflection, fi <H- —fi. This 
situation is realized if the quantity in question has a nonanalyticity along the #-axis, induced by 
the RW transition. Examples displaying this behavior are the quark number density n q and the 
imaginary part Im \P of the modified Polyakov loop IllOn . The C symmetry is explicitly broken 
whenever 9^0 and 9 ^ (2n — l)n/3, but close to the RW endpoint this symmetry can still be 
regarded as approximate. 

This discussion clarifies that the RW transition lines are first order, but it does not identify 
the order of the transition right at the RW endpoint: both first- and second-order transitions are 
possible at that point. In our nonlocal PNJL model, the first-order transition will indeed be shown 
to proceed on to the RW endpoint. In this case, the RW endpoint can become a triple point at 
which three first-order-transition lines meet. However, as 9 moves away from (2n — l)ir/3, the 
explicit C-symmetry breaking takes over and the first-order phase transition turns into a crossover. 

It has been suggested that the deconfinement crossover at fj, — can be considered as the 
remnant of the C-symmetry breaking (see Ref. [10] for details). As already mentioned, one can 
understand from this viewpoint that the transition behavior at finite n\ is connected to the chiral 
and deconfinement transitions at finite //r. Investigations of the [i\ region have indeed attracted 
considerable attention lately. Much recent progress has been achieved by analyzing the RW end- 

nn 

point in LQCD simulations with two- and three-flavor quarks as reported in Ref. [49, 50] and 
Ref. [26]. (We mention, in passing, that the RW endpoint is also investigated within the frame of 
gauge/string duality 115 ill ). 
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F. Nonderivative vector couplings 



The standard NJL and PNJL approaches usually work with a chirally symmetric combination 
of scalar and pseudoscalar quark couplings as a minimal setup. On the other hand, it is well 
known that additional vector and axial- vector couplings also emerge as parts of the chirally sym- 
metric effective interactions between quark quasiparticles. Here we shall focus on the role of 
isoscalar-vector interactions and their effects at imaginary chemical potentials. In the mean-field 
approximation this isoscalar-vector interaction is directly related to nonzero baryon density; it ap- 
pears independently of the derivative coupling already mentioned [the one that generates the quark 
wave-function-renormalization factor Z(p)]. As elaborated in detail in the Appendix, such a non- 
derivative vector interaction arises naturally by Fierz transformation from the interaction between 
quark color currents that has in turn its origin directly in QCD. 

The actual form of the additional nonlocal vector interaction used here is 

5C int = -G v f(x) 3fl (x), (11.24) 
f{x) = J d 4 z C{z) q(x + z/2) 7 M q(x - z/2) . (11.25) 

Introducing a mean field u associated with the Euclidean time component (density) of the current 
(III.25I) after bosonization, the primary effect of this interaction is to shift the chemical potential in 
Eqs. (IILT41) and flLffl) as 

li — > fi — u , (11.26) 
while the thermodynamic potential receives an extra contribution 

n -+ n - ^- . (ii.27) 

The coupling strength G v of the vector interaction is conveniently expressed in terms of the ratio 
Gy/G, relative to the scalar-pseudoscalar coupling G that controls spontaneous chiral symmetry 
breaking. In the actual calculations this ratio will be allowed to vary within 0.25 < G v /G < 0.5, 
corresponding to limiting axial U(1)a anomaly scenarios as discussed in detail in the Appendix. 

Note that the nonlocality distributions C(p 2 ) and F(p 2 ) are introduced at the Lagrangian level; 
hence they are to be calculated with chemical potentials // not shifted by the vector mean field to. 
Note, furthermore, that the thermodynamical potential must be real, so extensions to imaginary 



chemical potential imply imaginary to at the same time llllh 
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HI. RESULTS 



As we proceed to the actual calculations, the following versions of the nonlocal PNJL model 
are used in comparison: 

• Set I: without a Z factor [i.e., Z(p 2 ) = 1] and setting G v = 0. 

• Set ll:Z=l but G v ^ 0. 

• Set III: Z(p 2 ) included but G v = 0. 

• Set IV: Z(p 2 ) and G v ^ both included. 

This section first presents and discusses numerical results at imaginary chemical potential us- 
ing the nonlocal PNJL model in its versions with sets I and II, i.e., setting the wave-function- 
renormalization factor Z = 1 and comparing scenarios without and with vector interactions. In 
the second part, using sets III and IV, it will be demonstrated that additional effects induced by the 
Z{jp 2 ) factor are indeed individually small but can produce nonnegligible effects in the T-9 phase 
diagram when combined with vector couplings. 

Figure [2] shows the temperature dependence of the scalar mean field a, and of the real part 
of the modified Polyakov loop (111.21 1) , at different values of the imaginary chemical potential: 
= fr/T = 0, vr/6 and vr/3. The transition temperature at 9 = is T c = 208 MeV. At 9 = and 
7r/6, the a field and Re \P show the typical features of a crossover. At 9 = n/3 these quantities 
have already developed discontinuities, displaying gaps in both a and Re at T = 240 MeV. This 
means that the RW endpoint becomes a first-order transition point associated with spontaneous C- 
symmetry breaking. This feature repeats itself when 9 is shifted to 9 + 2irn/3, reflecting the RW 
periodicity. Our results at this point are consistent with LQCD data using finite-size scaling and 
two-flavor staggered quarks 

The transition temperatures for the chiral and deconfinement crossovers, T c x and T c d , have been 
determined at the maximum slopes of the order parameters. It is found that these pseudocritical 
temperatures coincide both at /i = and at finite jjj. This implies that two crossover lines are 
connected at the RW endpoint. Both lines can be considered as remnants of spontaneous C- 
symmetry breaking at 9 = ir/3. This interesting property is specific to the nonlocal PNJL model. 
In local versions of this model, with an artificial cutoff in momentum space, T c x and T c d may differ. 
In this case the chiral crossover line does not connect with the first-order RW endpoint, and the 
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Fig. 2. The T dependence of the scalar mean field a and of the real part Re \P of the modified Polyakov 
loop. The dotted, dashed and solid lines denote PNJL results at 8 = 0, 7r/6 and 7r/3, respectively, using 
input set I. 

chiral transition cannot be associated with a remnant of spontaneous C-symmetry breaking. Such 
a feature is also observed in the chiral limit at finite fi\. In this case the chiral transition is second 
order and the transition line terminates at a point different from the RW endpoint; see for example 



Fig. 1 in Ref. Ml 311. 



Next we examine the region in the vicinity of the RW endpoint at 9 = 7r/3. Figure [3] shows 
the T dependence of the imaginary part Im \P of the modified Polyakov loop in the nonlocal PNJL 
model. The dotted, dashed, and solid lines show results obtained at selected values 9 = 7r/3, 
147r/48, and 157r/48, respectively, of the imaginary chemical potential divided by temperature. 
One observes that Im \P, the order parameter of spontaneous C-symmetry breaking, develops a 
gap near 9 = n/3. Other quantities such as a and Re \P also display gaps at the same position, 



52D. However, 



reflecting the coexistence condition for nonanalyticities of order parameters f | 1 3L 
this first-order behavior starts being suppressed as soon as 9 deviates from 9 = ir/3. This be- 
havior is associated with explicit C-symmetry breaking becoming strong below the RW endpoint 
temperature Te- 

Figures H] (a), (b) and (c) display the 9 dependences of the order parameters a, Re W and Im \P, 
again calculated using the nonlocal PNJL model with set I. The dashed, dotted and solid lines 
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Fig. 3. The T dependence of Im ^ in the vicinity of the RW endpoint. The dotted, dashed and solid lines 
are nonlocal PNJL results (using set I) at 9 = it/ 3 and its close neighborhood. 



are the results at T = 220, 230 and 250 MeV, respectively. At T = 220 and 230 MeV, these 
quantities show smooth behavior around 9 = n/3. At T = 250 MeV, on the other hand, the 
#-even functions a and Re \P have sharp cusps at 9 = ir/3. At the same temperature, the 9-odd 
quantity Im & develops a characteristic gap at 9 = n/3. This is the RW transition mentioned 
earlier. These results confirm that our nonlocal PNJL model is indeed capable of reproducing the 
RW periodicity as well as the RW transition in a way consistent with QCD. 

It is important to note that, in the nonlocal PNJL model, the RW properties just described 
emerge from the consistent insertion of the Polyakov-loop field in the distributions, C(p 2 = co 2 i + 
p 2 ) and J-'ip 2 = to 2 i + p 2 ), with u nji given by Eq. (III. 141) . Had we omitted the combinations of 
gluonic background fields A\ and A\ in the shifted Matsubara frequencies cu n ^ , the extended Z 3 
symmetry would have been explicitly broken, with the consequence of losing the RW periodicity 
and maintaining only the 2n periodicity in 9 = jXi/T. The dependence of the chiral four-quark 



3411 . is thus a 



coupling on the A± 8 fields generating the Polyakov loop, as discussed in Refs. H3 1L 
natural consequence of the RW periodicity requirement. The nonlocal PNJL approach makes this 
important relationship explicit. 

Figure [5] shows the phase diagram at finite 9 in the nonlocal PNJL model. The solid lines 
represent first-order phase transitions while the dotted lines describe crossovers. The symbols are 
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Fig. 4. The 9 dependences of the scalar mean field a [subfigure (a)], and of the modified Polyakov-loop, 
Re & and Im & (subfigures (b) and (c)). The dashed, dotted and solid lines are results of nonlocal PNJL 
calculations (set I) at temperatures T = 220, 230 and 250 MeV, respectively. 



LQCD data taken from Refs. [22, 



25J] [here we translate p(a)/p(0) of the LQCD data to T/T c 



using the two-loop perturbative solution to the renormalization-group equation relating the lattice 
spacing a and the lattice gauge coupling (3 (a)]. The RW endpoint becomes a triple point at which 
three first-order transition lines merge. In this nonlocal PNJL model the transition lines determined 
by a and Im \P coincide in the entire 9 region. This confirms once again that the two crossover 
lines to the left and right of 9 = n/3 are remnants of spontaneous C-symmetry breaking in this 
model. 
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Fig. 5. Phase diagram of the nonlocal PNJL model (set I) in the T-9 plane, in comparison with LQCD 
results. The dotted and solid lines represent crossover and first-order transitions, respectively. Open squares 



and closed circles are LQCD data taken from Refs. [22] and Ref. I25I1 . The temperature is given in units 



of the pseudo-critical temperature T c at \i = 0, while the imaginary chemical potential is expressed as 
9 = m/T in units of 7r/3. 



Next, we turn to set II and investigate the effects of the additional vector interaction, Eqs. (III. 24!) 
and (III.25I) . as a function of the vector coupling strength G v . Figure [6] shows the G v dependence 
of the order parameter lm\l/ at 9 = n/3. The primary effect of the additional vector-type four- 
quark interaction is evidently a reduction of the temperature at which the triple point appears. This 
shift is welcome recalling the comparison with LQCD data in Fig. |5] At this stage a coupling 
G v /G ~ 0.6, slightly above the preferred range G v /G < 0.5 (see the Appendix), would give good 
agreement with lattice results. 

In the next step we examine the influence of the quasiparticle Z factor using the set III version 
of the nonlocal PNJL model. Figure|7]shows the T dependence of the order parameters at 9 = n/3 
and 9 = 0, with and without wave function renormalization. Evidently, the Z factor does not affect 
the order parameters in any significant way over the entire range of T and fi. 

Finally, moving to set IV the combined effects of the additional vector interaction and of the 
Z factor are studied. Figure [8] shows the dependence of Im \P on the vector coupling strength G v . 
The additional role of the Z factor is now to enhance the vector coupling effect in the downward 
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Fig. 6. The T dependence of Im & in the nonlocal PNJL + G v model at 6 = it/ 3. The dotted, dashed and 
solid lines are results obtained with G v /G = 0, 0.4 and 0.6, respectively. 
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Fig. 7. The T dependence of the order parameters a and Re ^ at 9 = and 6 = tt/3 in the nonlocal PNJL 
model (set III) with inclusion of the quark quasiparticle wave-function-renormalization factor Z(p 2 ) (solid 
curves) as compared to the case with Z = 1 (dotted and dashed lines). 
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Fig. 8. The G v dependence of Im & in the nonlocal PNJL model (set IV). The dashed and solid lines denote 
the results at 9 = tt/3 with G v = and G v = 0.4 G, respectively. 



shift of the temperature at which the triple point appears. In fact, good agreement with LQCD 
results is now achieved using G v /G = 0.4, well within the natural range of G v implied by the 
Fierz-transformed color-current-current interaction between quarks. The resulting phase diagram 
is displayed in Fig. [9] along with LQCD data. A direct comparison must take into account the 
fact that these data have been generated with a quark mass m = 12 MeV, about 4 times the 
standard input current quark mass in the PNJL calculation. The pion mass corresponding to this 
heavier quark mass is thus twice as large as the physical one. We have studied the quark mass 
dependence in the vicinity of the triple point. The open circle in Fig. [9] is the PNJL result using 
m = 12 MeV. Evidently, the difference in T/T c between calculations with physical pion mass and 
with m w ~ 280 MeV is only marginal. 

The behavior of the thermodynamic potential Q as a function of the (imaginary) isoscalar- 
vector mean field, Im to, is also of some interest here in order to examine details of the phase 
transition. Figure[l0] shows Q for set IV, once again using G v /G = 0.4, at three neighboring 
temperatures, demonstrating the occurrence of a first-order transition near T = 217.5 MeV. 
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Fig. 9. Phase diagram of the nonlocal PNJL model with inclusion of the quark quasiparticle renormalization 
factor Z{p 2 ). Upper curves: without additional vector interaction (model set III). Lower curves: including 
vector coupling with G v /G = 0.4 (model set IV). The dotted and solid lines represent crossover and first- 
order phase transitions, respectively. Also shown for comparison is the result with G v /G = 0.4 but with no 



Z-factor effect (set II). Open squares and closed circles are LQCD data taken from Ref [22] and Ref. [25]. 
The temperature scale T c is the pseudo-critical temperature at \i = 0. The open circle at the triple point is 
the result of the set IV calculation with quark mass tuq = 12MeV. 



IV. SUMMARY 



In this paper we have investigated properties of two-flavor QCD at finite imaginary quark chem- 
ical potential /ij in order to establish constraints for the modeling of the QCD phase diagram. A 
recently extended version of the nonlocal PNJL model has been used for this task. Apart from 
the momentum-dependent quasiparticle mass M(p 2 ) of the quarks, and their gauge-covariant cou- 
pling to a Polyakov-loop background field, this model also incorporates the quark wave-function- 
renormalization factor, Z(p 2 ). Both M(p 2 ) and Z(p 2 ) are introduced in close correspondence with 
results from Dyson-Schwinger calculations and LQCD computations. We can draw the following 
conclusions: 

1 . The nonlocal PNJL model is capable of reproducing important generic properties of QCD 
at jii > 0: the Roberge- Weiss (RW) periodicity and the RW phase transition. It turns out that the 
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Fig. 10. Thermodynamical potential Q of the nonlocal PNJL model (set IV with G v = 0.4 G) as a function 
of the imaginary vector mean field Im oj. The lines are results of calculations at neighboring temperatures 
T = 216.5, T = 217.5 and 218.5 MeV, as indicated. The reference tt is defined at Im w = 0. 

RW endpoint located at odd multiples of 9 = H\/T = it/ 3 is first order. It is actually a triple point 
at which three first-order phase transition lines merge, signaling spontaneous breaking of charge- 
conjugation symmetry. This first-order transition region in the T-9 phase diagram is restricted to 
a small vicinity of the triple point, outside of which explicit C-symmetry breaking takes over and 
turns the transitions into continuous crossover lines extending symmetrically left and right of the 
triple-point area. 

2. The (pseudo-)critical temperatures for the chiral and deconfinement crossovers and first- 
order transitions turn out to coincide both at zero chemical potential and at fxj > 0. This coinci- 
dence is a characteristic feature of the nonlocal PNJL model. It is based on the consistent treatment 
of the Polyakov-loop dependence of the nonlocality distributions that govern the effective interac- 
tions between quark quasiparticles. 

3. Isoscalar- vector interactions between quarks, additional to the standard chiral combination 
of isovector-pseudoscalar and isoscalar-scalar interactions, play an important role in the T-fx phase 
diagram at imaginary /i. Such vector interactions are well known to emerge naturally from Fierz- 
transforming the basic color-current-current interaction of QCD. Their presence is found necessary 
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in order to reproduce the location of the RW endpoint on the T scale and the pattern of transition 
lines in the neighborhood of that point as observed in LQCD. 

4. The wave-function-renormalization factor Z(p 2 ), resulting from derivative couplings be- 
tween quarks, does not play an important role individually. However, in combination with the 
isoscalar- vector interaction just mentioned, it has a visible impact on the phase diagram. The Z 
factor and the vector coupling cooperate coherently in reproducing the phase diagram at \i\ > 
from LQCD. Including both effects the resulting vector coupling strength, relative to that of the 
scalar interaction that generates spontaneous chiral symmetry breaking, is found to be G v = 0.4 G. 
This value is remarkably consistent with earlier studies of meson properties and spectroscopy 
within NJL models including realistic constraints from the QCD axial anomaly. 



V. APPENDIX: 

COLOR-CURRENT-CURRENT CORRELATIONS AND ISOSCALAR- VECTOR INTERACTION 

Here we first sketch the derivation (see Ref. J^j ) of an effective color-current-current inter- 
action from two-flavor QCD and then recall its Fierz transform leading to NJL-type four-point 
couplings (see e.g. Refs. ID3I . I54I1 ). Consider the QCD Lagrangian 



C 



QCD 



C + Cyu + £ 



-qG 



C q = qix)^^ - m)q(x) , 



-- F ; u (x)Fr(x) 



£ym — 
£ q a = gq(x)j ll A li (x)q(x) 



(V.l) 
(V.2) 

(V.3) 

(V.4) 



with q(x) = (u(x), d(x)) T and N 1 = A a \ a /2. The QCD generating functional becomes 



QCD 



VqVq J VA*exp(i J d 4 x C QCD 



VqVq exp 



d 4 x(£ g + iW[j}) 



where 



iW[j] = In J VAl exp [ J d 4 x (£ YM + C qG 



involves the quark color currents 



Qi x ) "7^"7m q( x ) 



(V.5) 



(V.6) 



(V.7) 



21 



Expanding the generating functional in powers of this color current, the function iW can be ex- 
pressed as 



iW\j] = iW[0]-g J d A xW^ a {x)f a \ 

+ 7 / d ' x I d ' yj " {x) w !2*( x > y) ib{y) + °(f) , (v.8) 

where is related to the connected gluon n-point function without quark loops. Truncating 
this expansion at W^ 2 \ this becomes a model in which the local color gauge symmetry of QCD is 
reduced to a global SU(3) C symmetry [global color model (GCM)], with the reduced Lagrangian 

Cgcm(x) = q(x) (< 7 „0" - m) q(x) + g 2 ^(x) [ d 4 y W$ ab (x, y) %(y) . (V.9) 



This Lagrangian features the nonlocal color-current interaction that involves the full gluon prop- 
agator. The nonlocal PNJL approach basically keeps this structure. Assume now that the gluonic 
correlator determining is very short-ranged so that it can be localized with a leading term 
proportional to S 4 (x — y). This is the idea behind the classic NJL and local PNJL models. With 
this local approximation the interaction Lagrangian becomes 

8 

C iat (x) = G C J2 [?(*) t a l, q(x)] [q(x) *V ?(*)] 
o=l 

8 

= G C J2 € a (t a U (7/.)™ (t a ) 7 s (Y)tu S kl q l uS , (V.IO) 

a=l 

with the SU(3) C generators t a = X a /2. Here G c is an effective color coupling strength of dimension 
[length] 2 . The second part of Eq. (IV. 101) explicitly displays the combinations and contractions of 
quark color indices a, f3, . . . £ {1, 2, 3}, flavor (isospin) indices £ {1, 2} and Dirac spinor 

indices r,s, . . . £ {1, . . . , 4}. This interaction represents color-octet exchange between quarks. A 
Fierz transform translates this into various color-singlet channels and color-octet quark-antiquark 
channels with isospin singlet or triplet quantum numbers. 

Consider first the general SU(iV) Fierz transformation from channels (12 — > 34) to the ex- 
change channels (14 — > 32): 

(1 ,g, 1) 1234 = 1(1 ® l) 14i32 + 2{t a <g> t a ) 14 , 3 2 , (VI 1) 

it* <g> t a )l 2 ,34 = \ (l - J (1 ® 1)14,32 - ^it a <g> t a )l4,32 • (V.12) 
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The "1" symbolizes the N x N unit matrix while the matrices t a stand for the generators r a /2 
or \ a /2 of SU(2) or SU(3), respectively. Using this expression, the Fierz transform of the color 
degrees of freedom in Eq. (IV. 101) becomes 

1 



If J_ 

2 V 1 ~ iVf 



(V.13) 



where the first term denotes the color-singlet and the second one refers to the color-octet chan- 
nels. In the present study we consider only the color-singlet part and absorb the coefficient 4/9 
multiplying (1 (g> 1) in a redefined coupling constant G c . 

Next, consider the Fierz transformation in Dirac space. The well-known result is 



(7m ® Y 



lrs,tu 



-(1(8)1 + 275 <g> ry 5 )r«,ts + ^(7^ ® 7 M + 7/^75 ® 7 M 75)r«,fe 



(V.14) 



Finally, the Fierz transformation in flavor-SU(2) space, following Eq. (|V.l II) . gives in terms of the 
isospin Pauli matrices t 1 : 

1 1 3 

(1 <g> l)y )fc , = -(1 ® l) <l)fcj - + - ^(r* <g) r'^jy . (V.15) 



i=l 



The color-singlet-interaction Lagrangian derived from the color-current-current interaction be- 
comes 

G 
2~ 

' ' ' ' ' 11 (V.16) 



C 



(o) 



(<?<?) 2 + {qilbqf + (<?-rg) 2 + (qij 5 rq) 2 

- ^{(<?7 M 9) 2 + (g7 M 7 5 g) 2 + (q^rq) 2 + (g7 M 7 5 Tg) 2 

This interaction still has a chiral U(2) R x U(2) L symmetry. An axial U(l)-breaking term must 
be added in order to account for the axial anomaly in QCD and reduce the symmetry to chiral 
SU(2) R x SU(2) L , times U(l)y for conserved baryon number. The U(1)a anomaly effect is 
introduced by the Kobayashi-Maskawa-'t Hooft determinant interaction as 



^anomaly = G A det[<? (1 + 75) q] + h.C. 

Ga 



(qq) +(qil5rq) -{qfq) -(qirf 5 qy 



(V.17) 



where det acts in flavor space. Setting Ga = G leads to the usual NJL-type Lagrangian (see 



Refs. [153 



5611 for reviews of the NJL model). In this case the total interaction becomes 



C 



int. 



£ (0) + C 



anomaly 



G 



(qq) +(qilsTq) d 



~ i") (qi^qf + (qi^qf + {q-ffqf + (qj^rqf 



(V.18) 
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This is the situation with "maximal" U(1)a symmetry breaking in which the anomaly term is 
strong and accounts for half of the scalar-isoscalar and pseudoscalar-isovector interactions, while 
the scalar-isovector and pseudoscalar-isoscalar terms are completely eliminated. In this case the 
vector coupling strength is G v = 0.25 G, denoting the vector interaction term as — G v (q^^q) 2 . 
The instanton liquid model IptI l^l produces such a scenario. However, the strength of the 
anomaly-driven part of the interaction is subject to ambiguity. In Refs. [0,0] the percentage 
for the anomaly term in the total scalar-isoscalar and pseudoscalar-isovector coupling constant 
was analyzed using the three-flavor NJL model. According to these evaluations, the anomaly 



term accounts for about one fifth (more precisely: 



interaction with parameters as given in Refs. B5 



61, 



6 %-21 %) of the total effective four-quark 



6211 . In this case the resulting ratio of vector- 



to-scalar couplings is G v /G = 0.4. In practice the acceptable range for this ratio is taken to be 
0.25 < G v /G < 0.5 with a preference for values around 0.4. 
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